Cosmological constraints from CMB distortion 
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We examine bounds on adiabatic and isocurvature density fluctuations from /i-type spectral dis- 
tortions of the cosmic microwave background (CMB). Studies of such distortion are complementary 
to CMB measurements of the spectral index and its running, and will help to constrain these pa- 
rameters on significantly smaller scales. We show that a detection on the order of fj, ~ 10 -7 would 
strongly be at odds with the standard cosmological model of a nearly scale-invariant spectrum of adi- 
abatic perturbations. Further, we find that given the current CMB constraints on the isocurvature 
mode amplitude, a nearly scale-invariant isocurvature mode (common in many curvaton models) 
cannot produce significant ^-distortion. Finally, we show that future experiments will strongly con- 
strain the amplitude of the isocurvature modes with a highly blue spectrum as predicted by certain 
axion models. 
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I. INTRODUCTION 

The scale-invariant formulation of the power spectrum 
of primordial fluctuations is highly successful [T] . The ob- 
servational results of both cosmic microwave background 
(CMB) and large scale structure (LSS) studies are mod- 
eled superbly by a scale-invariant spectrum. However, 
observational data for the power spectrum on small scales 
(smaller than 0.1 Mpc) is lacking, prompting discussions 
of a possible deviation from scale invariance on small 
scales. For example, many inflationary models predict 
a running spectral index, which may develop a large de- 
viation from a scale invariant spectrum at large CMB 
multipoles I (for a survey of theoretical motivations and 
observational prospects for a running index, see [2]). 

Isocurvature fluctuation modes can also produce de- 
viations from a scale invariant spectrum. For example, 
isocurvature fluctuations are generated in the axion dark 
matter model [3] and the curvaton scenario [4]. If the 
spectrum is blue, the isocurvature modes contribute ad- 
ditional power to the density fluctuations on small scales, 
although the isocurvature mode remains subdominant 
on large scales. Current observational constraints allow 
models with a mixture of scale-invariant adiabatic and 
isocurvature modes with a blue spectrum [5] , with isocur- 
vature modes constrained to contribute at most 10% on 
large scales (pure isocurvature models are already obser- 
vationally ruled out [5]). 

Currently, the running of the power spectrum index 
for adiabatic modes and the amplitude of isocurvature 
modes are constrained by CMB and large scale structure 
PQ. For the running power spectrum, the best fit value for 
WMAP and large scale structure is -0.022 ± 0.020. Fur- 
thermore, the amplitude ratio of the isocurvature mode 
to the adiabatic mode at k = 0.002Mpc~ 1 is 0.13 and 
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0.011 for the uncorrelated CDM isocurvature case and 
the correlated one, respectively. However these con- 
straints can be vastly improved by measurements of the 
fluctuations on smaller scales. 

The measurement of the CMB spectrum distortion 
from the perfect blackbody shape is a powerful tool for 
producing constraints on small scale density fluctuations 
before the recombination epoch. The energy of the den- 
sity fluctuation on small scales is released into CMB 
photons via Silk damping [7]. At a very high redshift 
(z > 10 6 ), the released energy is quickly thcrmalized by 
double-Compton and Compton scattering, and the spec- 
trum of the CMB photons maintains a blackbody shape. 
However, as the universe evolves, these scatterings be- 
come less effective and cannot establish thermalization 
of the released energy. As a result, released energy is im- 
printed on the CMB spectrum as CMB distortions from 
the blackbody shape [8] [9] . Such distortions may be clas- 
sified into two types, /x-distortion and y-distortion |10) . 
While ^(-distortion is produced from the double Compton 
scattering decoupling (z ~ 10 6 ) to the thermalization de- 
coupling by Compton scattering (z ~ 10 5 ), y-distortion 
is generated from thermalization decoupling by Compton 
scattering to the recombination epoch. In other words, 
by measuring these distortions, we can learn about the 
density fluctuations on scales smaller than the Silk damp- 
ing scale at these redshifts, up to k ~ 10 5 Mpc _1 . 

The current constraints on these distortions, obtained 
by COBE FIRAS, are < 9 x 10~ 5 and y < 1.5 x 
10~ 5 [11 . Proposals for future space missions such as 
PIXIE have the potential to provide dramatically tighter 
constraints on both types of distortion, with projected 
detection levels of ~ 5x 10~ 8 or y ~ 10~ 8 by PIXIE at 
the 5 a level [12] . It was recently pointed out that PIXIE 
can detect the CMB distortion due to the dissipation of 
the primordial fluctuations with n s — 0.96 at a 1 a level, 
because the distortions for n s = 0.96 are (i~8x 10~ 9 

[IMS]. 

In this paper, we investigate ^-distortions caused by 
the two types of fluctuations described above. One is the 
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primordial fluctuation with running spectral index and 
the other is an additional power-law fluctuation which 
is sub-dominant on large scales. One of the motivations 
for the latter type is the CDM or baryon isocurvature 
mode with a blue spectrum. An extremely blue spectrum 
(n-iso > 2) can be generated in the axion model suggested 
in 16J. Although the density fluctuations of CDM and 
baryons due to the isocurvature mode can survive Silk 
damping, the photon density fluctuations produced by 
the isocurvature mode suffer Silk damping |17j and this 
dissipation can produce /i-distortion. 



II. DISTORTIONS FROM CURVATURE 
PERTURBATIONS 

The energy injection in the early universe produces 
CMB spectral distortions. The evolution of the spectral 
distortions depends on the energy injection rate and the 
time scale of the thermalization process. In particular, in 
the regime between z ~ 10 6 and z ~ 10 5 , /z-distortion is 
created by the balance between the energy injection and 
the double Compton scattering. The time evolution of 
/i-distortion due to this energy injection is given by [3] 
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where tj^c is the time scale for double Compton scatter- 
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where Y p is the primordial helium mass fraction. The 
solution to this evolution is 
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The dissipation of acoustic waves in the photon-baryon 
plasma due to Silk damping injects energy into the CMB 
and causes the CMB spectral distortions. 

We can write the energy density perturbation of the 
acoustic waves as 



Q=^[c^ 7 (x) 2 ) + Mx) 2 )] 

The density and velocity perturbation of photons are 
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The power spectrum Pi(k), where the subscript i repre- 
sents 7 and v, is related to the primordial power spectrum 
by 



Pi{k) = A 2 (k)Pl 



(8) 



where the transfer function A,(fc) for modes well inside 
the horizon satisfies the relation 
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and 
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with the sound horizon r s given by 
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and kjj is the diffusion scale 
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In the above, R is the baryon energy density 
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n e is the free electron number density (before the recom- 
bination epoch, it is given by n e (z) = (nuo + 2n# e0 )(l + 
z) 3 = no(l + -z) 3 ), and a? is the Thomson scattering 
cross-section. During radiation domination the diffusion 
scale becomes 

k D =A D 1/2 (l + z)^ ; 

A D = - — - = 5.92 x 10 10 Mpc 2 . (14) 
135-ffo^r n e Q(TT 

Now, we can write the energy density in acoustic waves 
is 
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The primordial power spectrum can then be related to 
the power spectrum for the curvature in the co-moving 
gauge, C, via [H] 
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where R v is the neutrino energy density 
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The power spectrum of the curvature is parameterized as 
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with normalized amplitude — 2.4 x 10 9 , and pivot 
scale ho — .002Mpc _1 , the spectral index is ns, and the 
running of the spectral index is a s = dn s /d\nk. 

The energy release per unit redshift is given by 
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We then solve Eq.(|3]) while scanning over the a s 



n s 

plane using the values Y p = .24, VL^h 2 — .0224, along 
with the previously mentioned values for Aq, Ad, and 
toe- We have calculated /i for the injection interval Z\ = 
2 x 10 6 to Z2 = 5 x 10 , and the results are displayed in 
Fig.Q. 

The current WMAP7 bounds are given in Table 7 of 
PP. Roughly, for no tensor modes and no running (a s = 
0), n s = .968 ± .012; including running they find n s = 
1.008 ± .042 with a s = -.022 ± .02. Including tensor 
modes gives n s = 1.07 ± .06, a s = -.042 ± .024, and 
r < .49 (95%CL). The running we have plotted in Fig.Q 
is significantly smaller, showing the constraining ability 
of such a precise /x measurement. 

With regards to a measurement of the running, a s , 
Planck may be able to reach sensitivities on the order of 
\a s \ ~ .005, and we see that measurements of fi at the 
level of 10~ 8 , in tandem with complementary measure- 
ments of n s , would be at least competitive with such a 
measurement. 

As we have stated above, the measurement of fi may 
be comparable to or possibly supersede bounds on the 
spectral index or its running given by CMB observations 
such as WMAP and Planck. However, it should be em- 
phasized that the \x measurement would give informa- 
tion about such cosmological parameters on much smaller 
scales, and from effects occurring at much higher red- 
shifts. Thus, making such measurements is quite attrac- 
tive, and would allow us to reach previously unexplored 
distance and time scales for cosmological observations. 



III. DISTORTIONS FROM ISOCURVATURE 
PERTURBATIONS 
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FIG. 1. Values of /i generated by energy injection due to 
the damping of acoustic oscillations in the a 3 -n s plane for 
adiabatic modes with amplitude 2.4 x 10 -9 . 



where the subscript i is c for the CDM isocurvature and 
b for baryon isocurvature. We parameterized the power 
spectrum of S similarly to that of the curvature pertur- 
bations 
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where we have neglected the possibility of a running 
isocurvature spectral index. The value of the amplitude, 
A iso , is constrained by CMB data to be A iso < AAq. 

The power spectrum of the isocurvature perturbations 
for photon density fluctuations or velocity is given by 
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where A.u so which can be written, well inside the horizon, 
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where k eq is the wave number cosseponding to the Hubble 
horizon at the equality time. The energy release per unit 
redshift is given by 
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The existence of additional density fluctuations in the 
early universe is another possible heat source for the spec- 
tral distortion. Candidates of such additional density 
fluctuations are isocurvature perturbations. 

We define the isocurvature perturbations between pho- 
tons and CDM (or baryons) as 

S = Si-^5 7 , (21) 



W*) = (jf) e~ k2/kl - (26) 

Integration of this expression, as in Eq.([3|, will give 
the value of /x. In Fig[2j we show the results. For com- 
parison, we plot the analogous relevant quantities for the 
adiabatic perturbations in Fig(3] 

Compared to the perturbations on large scales, the per- 
turbations on small scales do not have enough time to 
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models [IB] which could possibly alleviate such a degen- 
eracy, there may be other sources of isocurvature modes, 
such as cosmic defects that can contribute to the degen- 
eracy. However, a future measurement of \l distortion, 
along with complementary measurements of (or stronger 
bounds on) isocurvature modes from upcoming experi- 
ments such as Planck, may break the degeneracy or more 
strongly constrain the scenarios which can produce the 
distortion. 



IV. DISCUSSION AND CONCLUSIONS 



FIG. 2. Values of /t generated by energy injection due to the 
damping of acoustic oscillations as a function of the isocur- 
vature spectral index, mi ao , for a range of values for the am- 
plitude of the isocurvature mode, Ai ao . The maximum value 
shown for Ai SO corresponds to the upper bound determined 
by WMAP of A iso < 2.4 x 1CT 10 . The shaded region has been 
ruled out by the COBE-FIRAS instrument. 
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FIG. 3. Values of \i generated by energy injection due to the 
damping of acoustic oscillations is shown as a function of a 
generic adiabatic spectral index, n a d, for a range of values 
for the amplitude of the adiabatic mode, A a d- This figure is 
shown for comparison with Fig[2]in order to see the disparate 
contributions to /i from adiabatic and isocurvature modes. 
The shaded region has been ruled out by the COBE-FIRAS 
instrument. 



grow until they re-enter the horizon. The density pertur- 
bations on small scales are suppressed as shown Eq. ( 24 ) 



and the energy of the acoustic waves is also small. As a 
result, //-distortion cannot give as strong of a constraint 
on the isocurvature perturbations, compared to the adi- 
abatic perturbations. 

With the introduction of possible isocurvature con- 
tributions to the //-type distortion, one sees that a de- 
generacy between the various sources of the distortion 
may arise. This is due to the fact that // is simply a 
scalar quantity and there may be various contributions 
to its value. Although there may be correlations between 
isocurvature and adiabatic modes in various inflationary 



The dissipation of acoustic waves due to Silk damp- 
ing creates CMB spectral distortions of the blackbody 
spectrum. In this paper, we have calculated //-type dis- 
tortions from energy injection into the CMB due to the 
dissipation of acoustic waves from Silk damping, focus- 
ing on both the nearly-scale invariant adiabatic pertur- 
bations with a running spectral index, and isocurvature 
perturbations. 

Fig. [I] is our main result for /i-distortion created by 
the adiabatic perturbations with a running spectral in- 
dex. Creation of this //-distortion depends on the inte- 
grated dissipation energy from the decoupling of the dou- 
ble Compton scatterings to the decoupling of the ther- 
malization by Compton scatterings. We found that the 
dominant contribution comes from the dissipation energy 
at k ~ 100 Mpc -1 in the case of the nearly-scale invari- 
ant adiabatic perturbations. The power spectra, which 
have the same amplitude at k ~ 100 Mpc -1 , generate the 
same /i-distortion, even though each spectral index and 
running index are different. As a result, a degeneracy 
arises between n s and a s as shown in Fig. [I] Although 
the current constraint on //-distortion by COBE FIRAS 
does not give a significant constraint on both n s and a s , 
the proposed PIXIE observer [T2] has the potential to 
give a much tighter constraint, regardless of detection 
or non-detection of //-distortion. For example, the de- 
tection of 5 x 10 -5 //-distortion by PIXIE will give the 
constraint, —0.004 < a s < for the scale invariant spec- 
trum (n s — 1). This constraint is tighter than what is 
obtained by WMAP by nearly an order of magnitude. 
PLANCK will be able to give a constraint comparable 
to PIXIE, although Planck's constraint is the result of 
the measurement of the fluctuations on the scales larger 
than k ~ IMpc -1 . Therefore, these measurements will 
be complementary. 

We also have evaluated //-distortion due to the isocur- 
vature perturbations. The isocurvature perturbations 
can grow outside the horizon. Therefore, the smaller the 
scale of the perturbations, the lower their amplitudes. As 
a result, //-distortions created by the isocurvature per- 
turbations are smaller compared to those generated by 
the adiabatic perturbations with the same amplitude and 
spectral index. We have found, considering the current 
constraint on the isocurvature amplitude from WMAP, 
a nearly scale-invariant isocurvature mode, which is mo- 
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tivated by the curvaton scenario, cannot produce signifi- 
cant ^-distortion. Concerning possible cross-correlations 
between adiabatic and isocurvature modes then, it is con- 
jectured that, even if a scale-invariant isocurvature mode 
is perfectly correlated with the adiabatic mode, the cor- 
rection to the /i-distortion compared to that arising solely 
from the adiabatic mode is very small. The /i-distortion 
measurements at the level proposed by PIXIE would then 
be inefficient in placing a constraint on the correlation 
between the scale-invariant adiabatic and isocurvature 
modes. However, PIXIE can place a strong constraint 
on the amplitude of the isocurvature mode with a very 
blue spectrum (n < 2.5), which can be generated, for 
example, in the axion model |16j . 

Finally, we comment on other potential sources of /i- 
distortion. The /i-distortion can be created by cooling 
of electrons, which was first computed in |13j . decaying 
particles [19], dark matter annihilation [20], dissipation 
of the magnetic field energy [21] . Hawking radiation of 
primordial black holes [H] , and cosmic strings [53]. In 
particular, /i-distortion due to the cooling by electrons is 
predicted in the standard cosmology and enters with a 
negative sign. Therefore, this distortion can cancel some 
or all of the positive /i-distortion due to the dissipation of 
acoustic waves. The predicted /i-distortion is —2 x 10~ 9 . 

In the present work we have used approximate analytic 
forms for the transfer functions, and therefore a more 
precise analysis is desirable. Although we expect the size 



of the contributions from such an analysis to be negligible 
at the sensitivity of PIXIE, we require a more systematic 
analysis of /i-distortion as PIXIE has the potential to give 
a precise constraint on the primordial perturbations. To 
address these issues, work is in progress to calculate the 
/i-distortion more precisely using numerical methods. 

As was point out in |14| , even a non-observation at 
the level of [i ~ 0(1O~ 9 ) would be a sensitive cosmo- 
logical probe. Furthermore, it should also be noted that 
if a /i-distortion were observed on the order of ~ 10~ 7 
or greater, the conventional cosmological picture would 
struggle to accommodate such a large value (see Fig[I]). 
Thus, such an observation would lead to the tantalizing 
possibility of new physic^] Taken together these points 
underscore the exciting insights that may emerge from 
future CMB blackbody measurements. 
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